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Abstract. New method of introducing vortex lines of the electromagnetic field 
is outlined. The vortex lines arise when a complex Riemann-Silberstein vector 
(E + iB)/\/2 is multiplied by a complex scalar function <j>. Such a multiplication 
may lead to new solutions of the Maxwell equations only when the electromagnetic 
field is null, i.e. when both relativistic invariants vanish. In general, zeroes of the 
<j> function give rise to electromagnetic vortices. The description of these vortices 
benefits from the ideas of Penrose, Robinson and Trautman developed in general 
relativity. 
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1. Introduction 

Complex scalar fields in a three-dimensional space carry generically a collection of 
vortex lines. These vortex lines are located at the zeroes of the field. Complex 
generic vector fields do not have such properties because the requirement that all field 
components simultaneously vanish leads to an overdetermined set of equations. There 
are at least two ways to overcome this difficulty and to introduce vortex lines also 
for vector fields. One may either select one relevant component of the field (as, for 
example, in the scalar theory of light) or one may build a single field entity from 
various vector components. The first approach has been used in most papers on phase 
singularities in wave fields (for a thorough review and list of references see pQ). In our 
previous publication [3] we have chosen the second method. We have studied phase 
singularities and vortex lines associated with the zeroes of the square F ■ F of the 
Riemann-SilbersteinJ (RS) vector (my units are chosen so that e = 1 and /i = 1) 



The vortex lines of the electromagnetic field defined in this manner are rather elusive 
objects — they lack clear signature and may be hard to observe. Even though the 
phase of the field has dislocations, the electromagnetic field does not vanish on these 
vortex lines, but it is just a null electromagnetic field — the two relativistic field 
invariants vanish on vortex lines. Nevertheless, the vortex lines built on F F may still 
have some role to play in singular optics as discussed in the papers by M. Berry [S] 
and G. Kaiser ,o| appearing in this issue. 

In the present paper I explore a different approach. It is based on the observation 
that when a null solution of the Maxwell equations is taken as the background field, an 
extra scalar multiplier may imprint on this solution a rich vortex structure. This time 
all components of the electromagnetic field will vanish on vortex lines. Of course, the 
vortex lines that are introduced in this manner do not have a generic character since 
they are found only in special cases — for null fields. However, the class of solutions 
of the Maxwell equations that carry these vortex lines is, in my opinion, sufficiently 
broad and interesting to make this approach relevant for singular optics. In addition, 
the present study touches upon some important concepts discovered before in the 
general relativistic context. 

2. Scalar prefactor as a carrier of phase singularities 

The starting point of this investigation is an observation that a complex scalar function 
(j)(r,t) multiplying the RS vector will control the zeroes of the electromagnetic field. 
Wherever </>(r, t) vanishes, the electric and magnetic field vectors vanish. Assuming 
that the background field F, as well as the product <f>F, satisfy the Maxwell equations 



f This name was introduced in a review paper |2J on the photon wave function, where one may also 
find information on the history of the Riemann-Silberstein vector. The first to make a practical use 
of the RS vector in the analysis of electromagnetic waves was Bateman g|. The RS vector offers a 
very convenient representation of the electromagnetic field, especially in the study of vortex lines, 
and I shall make an extensive use of it here. 




(1) 



id t F(r,t) = Vx F(r,t), V-F(r,t) = Q, 
id t (<t>(r,t)F(r,t)) = Vx ( ( f>(r,t)F(r,t)) , V-(<j>(r,t)F(r,t)) = 0, 



(2) 
(3) 
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we arrive at the following conditions on the prefactor field <f>(r,t) (to simplify the 
notation, from now on I omit the dependence of F on space-time coordinates) 

Fid t (j ) {r,t) = -FxV^r,t), (4) 

F-V0(r,t) = O. (5) 

This set of equations possesses nontrivial solutions only when the background field is 
null, i.e. F 2 = E 2 /2 — J3 2 /2 + \E B = 0. In order to prove this assertion, we can take 
the scalar product of both sides of Eq. (0} with the vector F. Since the scalar product 
on the right hand side is equal to zero, the left hand side must also vanish and this 
results in the following alternative: either dt<fi(r,t) or FF must vanish. In the first 
case, we end up with a trivial result: the scalar function (f)(r,t) must be a constant. 
It cannot depend on space and time variables since its gradient is at the same time 
parallel to F, as seen from Eq. and orthogonal to F, as seen from Eq. 0. In the 
second case we may have a nontrivial form of <p(r,t). Taking the scalar product of 
Eq. (j3J with the complex conjugate RS vector F* we obtain the condition 

iF* F dt<j>(r,t) = -F*-(F x V<f>(r,t)) = — (F* x F)-V<f>(r,t). (6) 

Finally, I obtain the following two partial differential equations for the function <fi(r, t) 

a^(r,t) + n-V^(r,t)=0, F- V</>(r, t) = 0, (7) 
where n is the normalized Poynting vector 

-iF* x F _ E x B 

n ~ F*-F ~ E 2 /2 + B 2 /2' ^' 

The solutions of Eqs. 10 will in general have zeroes and this will lead to vortex 
lines of the electromagnetic field riding atop the background solution F. On each 
vortex line the electromagnetic field vanishes. Near the vortex line the electric and 
magnetic field vectors followed around a closed contour rotate by 2itm, where m is the 
topological charge of the vortex. In contrast to nonrelativistic wave mechanics |H] , 
there is no interaction between vortex lines introduced in this manner. All vortices 
move independently, since the product of solutions of the first order partial differential 
equations J7J is again a solution. We may always add new vortices without changing 
the motion of the existing ones. In the following sections I shall give general solutions 
of Eqs. Q in two special cases: when the background field is a plane monochromatic 
wave and when it is the Robinson- Trautman field. 



3. Simple example 

Let us consider the solution of the Maxwell equations described by the following RS 
vector 

F(r, t) = (x + iy) exp(ifcz — not). (9) 

This solution describes the left-handed circularly polarized wave propagating in the z 
direction. It is a null field, since (x + iy) ■ (x + \y) = 0. In this case, Eqs. (J7| take on 
the form 

(d t + d z )cj)(r,t) = 0, {d x + id y )ct>(r,t)=0. (10) 

General solution of these equations is an arbitrary function of two variables: the real 
variable z — t and the complex variable x + iy. Thus, we arrive at the family of solutions 
of the Maxwell equations of the form 

F(r, t) =f{z-t,x + iy)(x + iy) exp(ifc(z - t)). (11) 
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These solutions may have a rich vortex structure. Assuming, for simplicity, a 
polynomial dependence on x + iy = w, we may write down the function / in the 
factorized form 

f{z — t,x + iy) = (w — ai(z - t))(w — a-i{z — t)) . . . (w — a n (z - t)), (12) 

where ak(z — t) are some functions of z — t. Each term in this product gives rise to an 
infinite vortex line whose shape is determined by the equation x + iy = ak(z-t). Each 
vortex line moves independently. In general, each vortex line flies with the speed of 
light in the z-direction. However, when the function ak(z — t) is just a constant, the 
corresponding vortex line is a stationary straight line in z-direction. A simple example 
of the 4> function of the type l|12fl is 

<t>(r,t) = (x + iy) m . (13) 

In this case the RS vector <f)F carries a straight vortex line along the z-axis with the 
topological charge m. 



4. TVautman-Robinson fields 



Null solutions of the Maxwell equations were extensively studied by Robinson [5] in 
connection with his work in general relativity. In particular, these solutions played 
an important role in the discovery of new solutions of Einstein field equations 10 
and also they appear in Penrose twistor theory A family of null solutions of the 
Maxwell equations that featured prominently in these studies was described in 12 - 
I shall call them the Robinson- Trautman (RT) fields§. They are usually expressed in 
a special coordinate system but I shall use here more intuitive Cartesian coordinates. 
In these coordinates the RS vectors representing the RT solutions can be written in 
the form 

F x = f{a,f}) /r l , (14a) 
t — m — z 

F y = if(a,(3) / 2 + 1 , (146) 
t — la — z 

F, = -f(a,0) - 2/3 , (14c) 
t — la — z 

where /(a,/3) is an arbitrary function of the following complex combinations of the 
Cartesian coordinates 

i • i x 2 +y 2 x-iy 
a = t — m + z : , p = ; . (15) 

t — la — z t — m — z 

The RT Maxwell fields are related to the simple plane wave field discussed in the 
previous section. Taking the / function in the form f(a,f3) — a -3 , we obtain the 
following localized null solution of Maxwell equations 

F x — ((x - iy) 2 — (t — ia — z) 2 ) /d 3 , (16a) 

F y = i {{x - iy) 2 + (t-ia- z) 2 ) /d 3 , (166) 

F z = -2{x-iy)(t-ia- z)/d 3 , (16c) 

where d = ((t — ia) 2 — x 2 — y 2 — z 2 ). This field may be obtained from the plane 
wave solution of the previous section by the coordinate transformation (conformal 

§ Topical review has been recently published by Trautman 1131 
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reflection) x M — > x^ /x 2 accompanied by the shift in time by an imaginary constant 
— ia and evaluated in the limit of infinite wavelength, when k — > 0. The energy density 
of this field is equal to 

(a 2 + (t-z) 2 + x 2 +y 2 ) 2 , s 

F*-F = 2 i 5 . 17 

((a 2 - t 2 + x 2 + y 2 + z 2 ) 2 + 4a 2 < 2 ) 3 

From this formula we see that the electromagnetic field ll(jalll6ell can never vanish. 
Note, that the presence of ia eliminates the singularity on the light cone and makes 
the energy of the field finite, 

d3rF *- F ^^- (18) 

There are some simple solutions with a single vortex line built on this background 
function. The simplest ones are obtained by taking the scalar multiplier function </> 
equal to f3 m . These solutions have a stationary vortex line along the ^-direction with 
the topological charge equal to —m. A bit more complex vortex line is obtained when 
<p = a + Aa(3. The vanishing of (f> leads to the following two equations for the vortex 
coordinates 

3 

2y = t, (x + 2a) 2 + z 2 = 3a 2 + -t 2 . (19) 

Thus, at each time the vortex line forms a circle lying on the uniformly moving y = t/2 
plane. The radius of the circle is shrinking until t = and then it starts expanding. 
Possibilities for constructing solutions with more intricate vortex lines are unlimited. 
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